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Abstract 

We consider a system of particles with arms that are activated randomly to grab other 
particles as a toy model for polymerization. We assume that the following two rules are 
fulfilled: once a particle has been grabbed then it cannot be grabbed again, and an arm 
cannot grab a particle that belongs to its own cluster. We are interested in the shape of 
| a typical polymer in the situation when the initial number of monomers is large and the 

numbers of arms of monomers are given by i.i.d. random variables. Our main result is a 
£N) \ limit theorem for the empirical distribution of polymers, where limit is expressed in terms 

of a Galton- Watson tree. 
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1 Introduction 

> : 

^ . We consider a system of interacting particles which might be used as a toy model for the 
formation of certain polymers. Roughly speaking, each particle possesses a certain number of 
arms that are activated successively at random to grab other particles. When an arm becomes 
active, it creates an oriented edge from the particle with the activated arm to the particle 
which has been grabbed, so the system describes the evolution of clusters of vertices with some 
inactive arms which are connected by oriented edges. As time passes, the number of inactive 
arms diminishes whereas the connectivity in the system increases by the creation of new edges. 
We impose the rules that once a particle has been grabbed, then it cannot be grabbed again, and 
that each time an arm is activated, it grabs a particle which does not belong to its own cluster. 
The latter requirement impedes the formation of cycles; thus the structure that connects each 
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cluster of particles is given by a tree and the terminal state obtained when all arms have been 
activated can be described ClS db random forest. 

This grabbing particle system bears some resemblance with the so-called random configu- 
ration model, which has been introduced independently by Bollobas [2] and Wormald [9]; see 
also Chapter 3 in Durrett |3J for further references. The latter is a natural stochastic algorithm 
which aims at producing a random graph on a set of vertices with pre-described degrees (in 
general the resulting graph is not simple and may contain loops and multiple edges). Typi- 
cally, a certain number of stubs is appended to each vertex, and one joins pairwise the stubs 
uniformly at random to create edges between vertices. Similarities and differences between the 
two models will be discussed at the end of this work. 

We are mainly interested in the case when the number of particles is large, and more precisely, 
in the situation when the empirical distribution of the number of arms converges to some 
probability measure \i on Z + as the number of particles goes to infinity. Thinking of these 
particles as monomers, our goal is to get information about the statistics of a typical polymer 
structure when the polymerization process is completed. In this direction, we shall consider the 
slightly different setting where the sequence of the number of arms is given by i.i.d. random 
variables with law /i, which is simpler to analyze. Loosely speaking, we will show that when the 
number of particles goes to infinity, a typical tree (i.e. a tree chosen uniformly at random in the 
terminal forest) tends to be distributed as a Galton- Watson tree with reproduction law given 
by the distribution of the number of arms of a typical particle. This result is perhaps quite 
intuitive; however providing a rigorous argument is not so straightforward, and this motivates 
the present work. 

The remainder of this paper is organized as follows. The grabbing particle system is presented 
formally in the next section, in the situation when the number of particles is finite and the 
number of arms of each particle is deterministic. In Section 3, we consider the case when 
the number of arms is random and given by a sequence of i.i.d. variables, and we make the 
connexion with Galton- Watson processes. More precisely, we shall show that conditionally on 
the total number of arms, the distribution of the terminal state is that of a Galton- Watson 
forest with a given number of trees, conditioned by its total size. Finally Section 4 is devoted 
to the asymptotic study of the empirical measure of trees in the terminal state when the number 
of particles tends to infinity. 

2 The grabbing particle system 

Consider n > 2 particles which are labeled by 1, . . . , n. Each particle, say i, consists in a 
vertex to which a certain number Xi G Z + of arms are attached. Arms can be either active 
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or inactive. An active arm attached to % is an oriented edge % — > j linking % to another vertex 
j 7^ i, whereas an inactive arm is an incomplete oriented edge % — > where only the left-end is 
specified. Initially, every arm is inactive, and once an arm has been activated, it remains active 
forever. We assume that the total number of arms is less than the number of particles, i.e. 

Xi + ■ ■ ■ + x n = n — k for some integer 1 < k < n. (2.1) 

We then enumerate these n — k arms uniformly at random, which specifies the order at which 
they are activated. 

Specifically, suppose that the first arm to be activated belongs to the particle i. Then at 
time 1, this arm grabs a particle j ^ % chosen uniformly at random amongst the other n — 1 
particles. The system at time 1 thus consists in n — 2 isolated particles with inactive arms, and 
a cluster of two particles, i and j, connected by an oriented edge i — > j which is labeled by 1. 
The particle i has still Xi — 1 inactive arms, and the particle j has kept its Xj inactive arms; 
inactive arms are labeled by 2, . . . , n — k. We iterate in an obvious independent manner, with 
the rules that every particle can be grabbed by at most one arm, and that each time an arm is 
activated, it grabs a particle from a different cluster. 

It is convenient to depict the evolution of the system in the upper half-plane by introducing 
some further orderings which of course do not alter the dynamics. Specifically, imagine that 
initially, particles lie on the horizontal axis, such that the initial order from the left to the right 
is uniformly random. We may then think of the X{ arms attached to the particle % as vectors 
making an angle £ir/ (xi + 1) for £ = 1, . . . , Xj with the (directed) horizontal axis. If, as above, 
the first arm to be activated belongs to the particle i and grabs particle j, then at time 1, the 
particle j (together with its arms) is shifted at the right-end of the first arm. See Figure 1 for 
this first step of the evolution of the system (with n = 7, k = 2, i = 6 and j = 5). 




Figure 1: planar representation of a grabbing particle system 
at times t = (left) and t = 1 (right) 



We then continue by an obvious iteration. Only particles lying on the horizontal axis can 
be grabbed, and each time such a particle is grabbed by the activation of an arm, it is shifted 
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together with its connected component to become the right-end of that arm. As time passes, 
we thus obtain a growing family of planar rooted trees, where vertices may bear some further 
inactive arms, and the roots are the vertices still lying on the horizontal axis, i.e. those which 
have not been grabbed so far. Edges and inactive arms are labeled by 1, . . . , n — k, the label 
of an edge being always smaller than that of an inactive arm, and the vertices are labeled 
by 1, ... ,n. The evolution reaches its final state at time n — k when all the arms have been 
activated; obviously we then have a forest of k planar rooted trees with n — k labeled edges on 
a set of n vertices, and the trees are ordered from the left to the right. See Figure 2 below. 
We stress that the terminal state retams a., the essenfafl information about evolution of the 
system, as it suffices to prune the last n — k — i edges to recover the clusters of particles in the 
system and their structures at time £. 



We shall now describe the distribution of the terminal state for the grabbing particle system. 
In this direction, it is convenient to introduce first the space $ nj fc of planar forests with k rooted 
trees on a set of n labeled vertices, with labeled edges. The k roots lie on the horizontal axis 
and are thus naturally ordered from the left to the right. The choice of a root in a tree induces 
an orientation of the edges in that tree, by deciding that each edge points away from the root. 
In this setting, recall the notation % — > j for i, j G {1, . . . , n} with i ^ j to denote an oriented 
edge from the vertex labeled by i to the vertex labeled by j. For any forest tp G § n ,k and 
j G {!,..., n}, we write 



for the outer-degree of the vertex labeled by i, i.e. the number of oriented edges which left-end 

1 However the exact ordering of particles lying on the horizontal axis is lost in general, but this is not relevant 
as far as empirical distributions of the clusters and their structures are concerned. 




Figure 2: terminal state as a planar forest 



di(ip) := Card{j G {!,..., n}\{i} : i — > j is an edge of (p} 
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is the i-ih vertex. 

Recall that xi, . . . , x n is a sequence of nonnegative integers such that (12. ip holds. We then 
denote by $(xi, . . . , x n ) C the subset of planar forests ip with k rooted trees, labeled edges 
and labeled vertices, such that the outer-degree of i is di(p>) = Xi for every i — 1, . . . , n. 

Lemma 1 For every integer 1 < k < n and every sequence of nonnegative integers x%, . . . , x n 
which fulfills (12.11) . the distribution of the terminal state of the grabbing particle system is given 
by the uniform law on $(xi, . . . , x n ). 

Proof: For i — 0, . . . , n — k, let $^(xi, . . . , x n ) denote the set of planar configurations which 
can be obtained from the grabbing particle system at time £, when initially the particle i 
possesses Xi inactive arms for i — 1, . . . , n. In other words, . . . , x n ) is the set of planar 

forests on a set of n labeled vertices, with n — i rooted trees lying on the horizontal axis, to 
which a total oin — k — i yet incomplete edges (i.e. with only a left-end) have been attached to 
certain vertices, in such a way that for every i — 1, . . . , n, the number of complete or incomplete 
edges with left-end i is Xi. Further, complete or incomplete edges are enumerated in such a 
way that complete edges are listed before the incomplete ones. In particular for i = n — k, in 
the notation above, there is the identity $ n _fc(xi, . . . , x n ) = $(xi, . . . , x n ). 

Plainly, under the assumptions made in the claim, the initial state has the uniform distri- 
bution in $0(^1? ■ ■ ■ 1 x n)- We claim that the activation of the first arm yields a random state 
in . . . , x n ) which has again the uniform distribution. Indeed, pick an arbitrary state 

ip G $i(xi, . . . ,x n ). There are exactly n states in $0(^1; • • • i x n) which may possibly yield ip 
by activating the first arm. More precisely, these are the states which are obtained from tp by 
pruning its unique complete edge at its right-end, and inserting this right-end (together with its 
inactive arms) on the horizontal axis relatively to the n — 1 roots of p. Further, starting from 
any of these n states, the probability of getting tp after activating their first arm is exactly the 
probability that this arm picks the right root amongst the n — 1 possible ones, that is l/(n — l). 
This establishes our claim. 

We then can iterate this argument, showing that the terminal state is uniformly distributed 
on $ n _fc(a;i, . . . , x n ) = . . . , x n ), as stated. □ 

3 Connexion with Galton- Watson forests 

We now turn our attention to the situation when the number of arms of particles is given by 
an i.i.d. sequence of random variables. Specifically, we consider a probability measure /i on Z + 
and assume that the number of arms of the i-th particle, which it is henceforth convenient to 
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denote by is random with law /i, and that £1, . . . , £ n are independent. We shall see that the 
terminal state of the system can then be described in terms of a Galton- Watson forest; in this 
direction, we first develop some material in this field. 

Consider a Galton- Watson process with reproduction law /i, so, roughly speaking, we have 
a population model with discrete non-overlapping generations, in which each individual gives 
birth at the next generation to a random number of children according to the law /i and 
independently of the other individuals. We shall work on the event that the total descent of 
each individual is finite a.s. By ranking the ancestors and the progeny of each individual, we 
can represent the process by a planar forest with rooted trees, where roots correspond to the 
ancestors at the initial generation, and each tree describes a genealogy, see e.g. Section 6.2 in 
Pitman [6]. 

We denote by F n ^ the finite set of all planar forests with k rooted trees on a set of n vertices 
(vertices and edges are not yet labeled). Recall that vertices and edges in a forest in F Utk 
can be listed by using, for instance, the depth-first search algorithm (cf. Figure 6.1 in [Sj on 
its page 126). In this setting, assigning labels to the n vertices (respectively, to the n — k 
edges) is equivalent to choosing a permutation on {1, . . . ,n} (respectively, a permutation on 
{l,...,n — k}). Thus, if for every I 6 N, E< stands for the group of permutations on {1, ... , £}, 
then there is the canonical identification as the product space 

^n,k Fyi.k ^ ^ ^n— k • 

We then introduce a product probability measure on F nk x E n x S n _ fc , denoted by GW^, that 
can thus be viewed as the law of a triplet of independent variables. The first random variable 
with values in F n fc is a Galton- Watson forest with reproduction law /i, having k trees and 
conditioned to have a total size n (assuming implicitly that this event has positive probability), 
and the second and third are simply uniformly distributed on E„ and respectively. We 

may thus think of GWj^ as the law on of a random state which is obtained by picking a 
planar forest with k rooted trees on a set of n vertices according to a Galton- Watson process 
with reproduction law /i, started with k ancestors and conditioned to have total size n, and 
independently enumerating the edges and the vertices uniformly at random. 

Theorem 1 Let £i, . . . be a sequence of i.i.d. variables with law /i. Fix integers 1 < k < n 
and work conditionally on £i + • • • + £ n = n — k, provided that the probability of this event is 
nonzero. Then the distribution of the terminal state of the grabbing particle system started with 
n particles and such that the initial number of inactive arms of the i-th particle is £j for every 
i G {1, . . . , n}, is the law GW^ on $ n ±- 

Proof: We first work with a Galton- Watson process with reproduction law \x started with 
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k ancestors; the notation P& will refer to the distribution of the process. Let (X{)i>o be a 
sequence of i.i.d. variables with law /i, introduce the downwards skip-free random walk St = 
Xi + - ■ - +Xi—£ and the first-passage time T := min{£ > 1 : Se = —k}. By listing individuals in 
the Galton- Watson process according to the depth-first search algorithm, it is well-known that 
Pfe can be identified with the distribution of the stopped sequence (Xi)x<i<T- More precisely, 
in this setting, Xi represents the number of children of the z-th individual listed by depth-first 
search. 

Next pick an arbitrary planar forest / G F nk with k rooted trees and n vertices (vertices 
and edges are not labeled). The depth- first search algorithm enables us to encode / by some 
sequence of integers {yi)\<i<n (more precisely, j/j is the outer-degree of the z-th vertex found by 
the depth- first search algorithm), and the probability that the genealogical forest induced by 
the k ancestors of the Galton- Watson process is given by / thus equals fi(yi) ■ ■ ■ fi(y n )- As a 
consequence, for every permutations a G E n and q G £ n -fc> we have 

with z = P fc (T = n). 

Then pick a planar forest (p with k rooted trees on a set of n labeled vertices, with labeled 
edges, which belongs to the subset $(xi, . . . , x n ). We can identify 

<p = (f, 

for a unique / G F nik , a G S n and q G S n _ fc . More precisely, / is encoded via depth- first search 
by the sequence i<*<n where X{ = y a U), and therefore 



n!(n — k)\z 



This shows that, provided that /i(xj) > for every i = 1, . . . ,n, the conditional probability 
GW^(- | $(xi, . . . , x n )) is simply the uniform law on $(xi, . . . , x n ). We deduce by comparison 
with Lemma [1] that for an arbitrary sequence of integers (xj)i<j<„ which fulfills (12. lft and 
fi(xi) > 0, the distribution of the terminal state in the grabbing particle system started with 
n particles and such that the number of arms of the z-th particle is Xi coincides with the 
conditional law GW<$(- | $(x x , . . . ,x n )). 

To conclude the proof, we just need to recall the following well-known consequence of the 
ballot theorem (see, for instance, Lemma 6.1 in Pitman [6]). Under the conditional law P/-(- | 
T = n) of the Galton- Watson process started with k ancestors and conditioned to have a total 
population with size n, the sequence of the number of children of individuals listed uniformly 
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at random is given by a sequence of n independent random variables with law n conditioned to 
add- up to n — k. Recall that the event that this sequence is given by (x±, . . . ,x n ) is precisely 
<&(xi, • • • , x n ). This completes the proof of the statement. □ 

4 Asymptotic behavior of the empirical measure 

Our goal here is to investigate the asymptotic behavior of the empirical measure of the trees 
resulting from the grabbing particle system started with n i.i.d. particles as n — > oo. More 
precisely, we shall work with the following setting. 

Let n be a critical or sub-critical probability measure on Z+, i.e. with mean 

oo 

m -=J2W) < !■ 

e=i 

We shall also always assume that /i(0) > 0. Let ^i, ... be a sequence of i.i.d. variables with law 
/x, and for each n > 2, set 

k{n) :=n -(£i + ••• + £„). 

We work conditionally on the event k(n) > 1 and consider the grabbing particle system started 
with n particles such that the initial number of inactive arms of the i-th particle is £j for every 
% G {1, . . . , n}. We write P n for the distribution of its terminal state in $ n := [J 1<k<n $ n ,k, and 
also denote by r±, . . . , Tk( n ) the sequence of planar rooted trees in this terminal state, ignoring 
labels on vertices and edges. That is to say that the sequence (ti, . . . , r fe („)) is the image of the 
terminal state by the canonical projection <& n ,k(n) F n ^ n y 

We also introduce the space of finite planar rooted trees with unlabeled edges and vertices, 

We denote by the probability measure on T which is induced by the genealogical tree of a 
single ancestor in a Galton- Watson process with reproduction law ji. We are now able to state 
our main limit theorem. 

Theorem 2 Assume that \x is a critical or sub-critical probability measure on Z + with /x(0) > 0. 
Then the proportion of tree structures identical to some fixed t e T in the terminal state of 
the grabbing particle system converges to T^(t) in probability under P n = P(- | k(n) > 1) as 
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n — > oo. T/iai is, equivalently, 

= (4.2) 



lim E n 



fe(n) 



We shall reduce the proof of (14.21) to establishing the following two asymptotic behaviors under 
the conditional probabilities P n . 

Lemma 2 The number of trees in the terminal state under the law P n converges in probability 
to oo, that is 

lim ~P n (k(n) > k) — 1 for every fceN. 

Proof: Note that 

(n) > fc) = P(6 + - • • + £n < n - k) 



tends to 1 when \i is subcritical (and then the statement is obvious), but may tend to when \i 
is critical. We therefore focus on the critical case. For simplicity, we assume that the support 
of \i is not contained into a strict subgroup of Z; the modification needed to treat the opposite 
case is straightforward. Then £ x + ' ' • + £n — n is an aperiodic, irreducible recurrent random 
walk on Z, and we know from a ratio limit theorem in Spitzer [7j on its page 49 that for every 
fixed t G Z 

P(6 + + = ~ P(£i + --- + £ n = n). 

It follows 

F(£ 1 + --- + £ n = n-£) = (P(e 1 + ... + e n <n)); 

consequently 

P(6 + • • • + £n < n - 1) ~ P(6 + ■ ■ ■ + < n) , (4.3) 
which yields our claim. □ 

Lemma 3 We have 

lim P n (n = t 1; r 2 = t 2 ,fc(n) > 3) = T^(ti)T^(t 2 ) . (4.4) 

n— >oo 

Proof: We denote by |t| the size (number of vertices) of a tree tel. Using Theorem [TJ we 
know that for each fixed n > 3 

Pn(Ti = ti, r 2 = t 2 , fe(n) > 3) 

n 

= c n J2nti + --- + Zn = n-k)GW i r l l l(T 1 =t l ,T 2 = t 2 ), (4.5) 

fc=3 
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with 

Cn = l/P(£ 1 + ... + £ n <n-l). 

Then recall the following celebrated identity of Dwass [4] : the first passage time process 

Tfc = inf{n > : £1 + ■ • • + £ n = n — k} 

has i.i.d. increments, which have the same distribution as the sequence of the sizes of trees in 
a Galton- Watson forest with reproduction law /i. In particular the step T\ is distributed as \r\ 
under T^. Further the formula of Kemperman (see, for instance, (6.3) in [6] on its page 122) 
states: 

k 

W>(T k = n) = -P(6 H h £„ = n - k) , for every 1 < k < n . 



By definition of the conditional law, we have thus for k > 3 and t := |ti| + |t 

GWi0t(Ti = ti,T2 = t 2 ) 



2 



T^(ti)T^(t 2 ) 



P(T fc = n) 

n{k - 2) P(^i H V £ n -t — n — k + 2 — t) 



k(n - t) P(6 + • • • + £ n = n - k) 
so combining with H4.5[) we arrive at 

Pn(n = ti, r 2 = t 2 , k(n) > 3) 
= Cn T^(t 1 )T M (t 2 )^- V -^-Hii + ■■■+ tn-t = n-k + 2-t). (4.6) 

fc=3 

In the subcritical case, we have c n —>■ 1 and also, by the weak law of large numbers, that the 
series in (14. 6 p converges to 1 as n — > oo. This proves (I4.4p in that case. 

We henceforth focus on the critical case. In this direction, from (14. 3 p and the identity 

oo 

P(6 + • • • + Cn+i < n + 1) = ^//(fc)P(ei + • • • + £„ < n + 1 - k) 

k=0 

we get by dominated convergence that 

P(£i + • • • + f„+i < n + 1) ~ P(6 + • • • + £ n < ra) . (4.7) 
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On the one hand, combining (j4.3j) . (14 .7p and the obvious upper-bound 
™ h — 2 

J2 — £— p (£i + • ■ ■ + = n- A; + 2-t)< P(& + ■ ■ • + f n _t < n - 1 - t) , 

fc=3 

we deduce from (14.61) that 

limsupP^n = tx,r 2 = t 2 ,k(n) > 3) < xW(ti)T M (t 2 ) . 

n^oo 

On the other hand, pick < t] < 1 and let k v be an integer such that (k — 2)/k > 1 — r] 
whenever k >k n . Combining (14.31) . (14.71) and the lower-bound 

n k-2 

— P (£i + • • • + Zn-t = n-k + 2-t)>{l- vMCi + ■■■ + £n-t <n-k v + 2-t), 

fc=3 

we deduce from (14.60 that 

liminf P n (n = ti,r 2 = t 2 , fc(n) > 3) > (1 - 7/)T^>(ti)T^(t 2 ) . 

Since 77 can be chosen arbitrarily small, this shows (I4.4j) . □ 

We have now all the ingredients to establish Theorem [2] using a standard argument of prop- 
agation of chaos (see e.g. Sznitman [8]). 

Proof: The obvious inequality 

Pn(ri = t 1 ,T 2 = t 2 ,k(n)>3) = P n (k(n)>3)<l= £ T^(t 1 )T^(t 2 ) 
ti,t 2 er ti,t 2 er 

combined with Scheffe's lemma show that the convergence (14.41) also holds in the space ^(TxT) 
of summable families indexed by the countable set T x T. It follows that 

lim P n (n=t) =T^(t), 

n—>oo 

and then by an argument based on exchangeability that 

/ , k(n) \ 
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Using again (14. 4p , exchangeability and Lemma [2J we get the asymptotic of the second moment 



Let us briefly discuss the super-critical case, i.e. when m > 1 (recall that we always assume 
that /i(0) > 0). Then F(k(n) > 1) tends to as n — > oo, so the conditioning on this event 
becomes singular. More precisely, it is readily seen that the distribution of k(n) under the 
conditional law P n = P(- | k(n) > 1) converges to some geometric law. In particular Lemma [2] 
fails and the situation differs crucially from that describes in Theorem [2j Nonetheless one can 
get an analogue of Theorem [2] in the supercritical case, provided that we condition on k(n) > cn 
for some < c < 1. Then standard arguments of large deviations can be incorporated to the 
proof of Theorem [2] and yield a similar limit theorem in which T^(t) has to be replaced by 
T^(t), where [L is the unique law in the exponential family of fx with mean 1 — c. Details are 
left to the interested reader (s). 

We now conclude this work by comparing the asymptotic behaviors for the present grabbing 
particle system and for the model of random configuration (see [H [5]), starting in both cases 
from a large number of particles such that the sequence of the numbers of arms attached to 
those particles is i.i.d. with a fixed law \i. Recall that in the latter, edges are formed by pairing 
arms uniformly at random, and that Molloy and Reed [5] have shown (under slightly different 
assumptions) that with high probability, the random configuration model does not possess a 
giant connected component if and only if Y^&Lo^ — 2)//(£) < 0. Further, in that case, if we 
select an arm uniformly at random and independently of the pairing process, then the structure 
of the cluster containing that arm can be described in terms of two independent Galton- Watson 
trees with reproduction law v(k) = (k + l)fi(k + l)/m; see Theorem 1 in pQ for details. So, 
roughly speaking, both models produce random graphs with the same degree sequences, and 
for both models the structure of a typical cluster can be described in terms of Galton- Watson 
trees. We stress that the reproduction laws are distinct (except in the Poisson case), which may 
seem paradoxical as degree sequences are closely related to the latter. However one should keep 
in mind that selecting a cluster by picking an arm uniformly at random differs from selecting a 
cluster uniformly at random, and more precisely induces a bias by the number of arms in the 
random cluster. This explains the relation between the two reproduction laws. 
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which establishes (14.21) . 



□ 
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have permitted to remove an irrelevant condition involving stable domains of attraction in our 
original version of Theorem 2. 
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